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Abstract 

It is shown that the quark-antiquark coalescence mechanism for 
pion production allows to explain the small pseudorapidity width of 
the balance function observed for central collisions of heavy ions, pro- 
vided effects of the finite acceptance region and of the transverse flow 
are taken into account. In contrast, the standard hadronic cluster 
model is not compatible with this data. 

1. Recently, measurements of balance functions (Ij in central collisions 
of heavy ions were reported by the STAR collaboration The striking 
feature observed in the data is the small width of the balance functions (in 
rapidity and in pseudorapidity), as compared to the expectations from the 
expanding thermally equilibrated quark-gluon plasma j2j . This indicates that 
hadronization occurs only at the very late stage of the development of the 
system 

One may then ask if the hadronization properties of the system produced 
in central collisions of heavy ions (as reflected in the balance function) are 
similar to those observed in nucleon-nucleon collisions. To investigate this 
problem we have evaluated the expected width of the balance function in 
pseudorapidity, using the pion cluster model which some time ago was suc- 
cessfully applied to nucleon-nucleon data 5J^. Corrections due to the finite 
acceptance region of and the effects of the transverse flow were included. 

*Address: Reymonta 4, 30-059 Krakow, Poland; c-mail:bialas@th. if. uj.edu. pi; 
^Note that in this model the quark-gluon structure of the system is entirely ignored. 
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Our estimates show that to obtain consistency with the data of [7 for 
central colhsions, the decay width of the pion cluster (in its rest frame) must 
be substantially narrower than that corresponding to isotropic decay. Thus 
one must conclude that the hadronization of the system produced in heavy 
ion collisions is rather different from that produced in hadronic collisions, 
where isotropic clusters can approximately account for the data 

Looking for a more adequate description, we considered the coalescence 
model j6j, which we generalized to include correlations. We thus assume 
that the hadronization proceeds in three steps. First, partons form neutral 
clusters. Subsequently, each cluster decays into some number of gluons and 
one^ qq pair (either uu or ddY- Quarks and antiquarks then recombine into 
positive, negative and neutral pions. The remaining gluons form new clusters 
and the process is continuing until all partons are transformed into hadrons. 

We show that this generalized coalescence model gives a good description 
of the data from j2], provided the decay of clusters into quarks and anti- 
quarks is isotropic, which seems a rather natural assumption. The obtained 
reduction of the width of the balance function (essential to account for the 
data) is a natural consequence of the coalescence process. It follows simply 
from the fact that the dispersion of the average of two independent random 
variables is smaller than the dispersion of each of them by factor \/2. 

It should be emphasized that we are discussing here only the angular 
distributions (expressed in terms of pseudorapidity) . The natural assumption 
of approximately isotropic and uncorrelated cluster decay is then sufficient to 
describe the width of the balance function. This is not the case for rapidity 
distribution where more detailed information on cluster decay is needed. 

Our conclusion is that the generalized coalescence model provides a natu- 
ral explanation of the very narrow width of the balance function observed in 
j2]. This result is a consequence of a very general fenomenon, characteristic 
for coalescence mechanism and thus, in our opinion, it does not depend on 
details of the specific structure of correlations between partons proposed in 
this note. Since, as we have seen, the model based solely on hadronic degrees 

^This conclusion is not surprizing since the measured rapidity width of the balance 
function in nucleon-nucleon collisions is about twice as large as that in central heavy 
ion collisions Our calculation shows that neither finite acceptance nor transverse flow 
effects can account for this difference. 

■^We consider just one pair for simplicity. This is not essential for the conclusions. 

^Throughout this paper by quarks and antiquarks we always mean -in the spirit of the 
coalescence model- the constituent quarks and antiquarks. 
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of freedom is not adequate, we feel that our result provides a rather strong 
argument in favour of the coalescence scenario. 

2. The balance functions can be expressed in terms of the single and 
double particle densities [SlIZI- Assuming, for simplicity, the (H — ) symmetry 
we have 

5(A2| Ai) = D{-, A2I+, Ai) - D{+, A2I+, Ai) (1) 

with 



Z}(+,A2|+,Ai) = D(-,A2|-,Ai) 



/a2 dV2 1 A, dv 



dr]\dri2 



■ dr]^ 



(2) 



/a2 dr]2 dr] 



d2 



n+_ 



D{-, A2I+, AO = D{+, A2I-, AO = ZT^"^^ (3) 

Jai dn+-^ 

where dn/drj and dn / driidri2 are the corresponding particle densities in pseu- 
dorapidity. 

The measurements of STAR require both particles to be in an interval 
—A < rji,rj2 < A (acceptance) while the difference of (pseudo) rapidities is 
kept fixed. This suggests a change of variables: 

Vi-V2 = 5; ivi + '72)/2 = z. (4) 

The integrations must be performed over dz with 6 being fixed. This implies 
that z must be kept inside the interval —A < z < A where A = A — \S\/2. 
Thus the balance function measured in [2^ is given by 



/_A dz [d'^n^ / dzd6 — d'^n^^ / dzdS] 

j:^^dr]+dn+/dr]^ 



3. Consider first a model in which pions are produced in neutral, isotropic 
clusters. It is well known that such models can account for the gross features 
of the nucleon-nucleon data [Hj. The distribution of clusters is denoted by 
Pc{y) where Y is the cluster pseudorapidity. 

To simplify the problem we assume that clusters decay into two charged 
particles and any number of neutrals. 
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The distribution in cluster decay is 

^ = M.-y)/(',.-y)/h--n (6) 

where h{z — Y) is responsible for correlations: /i = 1 if the decay products 
are uncorrelated. 

The single particle distribution in cluster decay p(r7+ — Y) is obtained 
by integration of (jH)) over rapidity of one particle. Both distributions are 
normalized to unity. 

The distribution of all particles is the convolution 



I dYp,iY)piv+ - Y) (7) 



with identical formula for negative particles. 

To evaluate the two-particle distributions one has to take into account 
that some pairs may come from different clusters and some others from one 
cluster. As is well-known (and can be easily confirmed by explicit calculation) 
the contribution from different clusters cancels in the balance function and 
thus only (H — ) pairs from one cluster do contribute. Their distribution is 
given by 



I dYp,{Y)h{z - F)/(77+ - - Y). (8) 



drj^drj^ 

Introducing (jHl) into (0) we have 

^ I-t dz J dYp,{Y)h{z - Y)f{v^ - Y)f{v_ - Y) 

j\dr^jdYp{r^-Y) ' ^ ^ 

4. To continue, we perform a simple exercise, assuming that all functions 
are Gaussians. We take^ 

Pc(r) = ^^exp(-yV^^); 
f{u) = J-expi-uyf); h{u) = yTTTW exp(-«V/) (10) 



^The normalization of the expression for h{u) guarantees the correct normahzation in 
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Using (fnH) one can evaluate the single particle distribution in cluster decay: 



This gives 



1 + /W 

P(^) = ^^exp(-r77a^); a = f \ . (11) 



< > -^exp(-5V2/^)-l=exp(-zVF2) (12) 
/v27r -T v^r 



with 



2[1 + /W] 



(13) 



This result introduced into ^ allows to express the integrals over z and y 
in terms of the error function. The result is 

which completes the calculation. 

This exercise shows that -in the limit of large acceptance- the width of 
the balance function is determined by the parameter / which describes the 
cluster decay. 

5. These results can be used to evaluate expectations from isotropically 
decaying clusters which were found roughly compatible with the data on 
hadron-hadron collisions [Sj- In this case the cluster decay distribution is 

^^^^ = 2(Jh^ 

giving the cluster decay width < |?7| >= log 2. This can be approximated 
by a Gaussian of the form ()11|) with a = log2A/7r ~ 1.23. Ignoring for the 
time beeing the effect of finite acceptance, we thus conclude from (fTT|) that 
the expected width of the balance function < > must be larger than 
a\plj ^ ^ .98 and thus by far exceeds the one measured in 0. 
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Finite acceptance [A = 1.3] of the STAR measurements j2] reduces the 
observed width of the balance function, as seen from (jl4|) . This is not suf- 
ficient, however, to bring the data in agreement with the model of isotropic 
pion clusters. The width < |5| > calculated from (fTT|) for isotropic clusters 
with f /g = (uncorrelated dcay) equals 0.67 at A = 3.5 (this value of A is 
roughly consistent with data for central collisions |HI), and does not change 
significantly when A varies around this value. Since, furthermore, < \6\ > 
increases with increasing f /g, there is no chance to meet the experimental 
value of 0.55. 

Another effect which may be responsible for the small width of the balance 
function is the transverse flow. Indeed, the clusters which are isotropic in 
their rest frame will not appear isotropic when moving with a transverse ve- 
locity. As shown in jOj , the distribution ()15|) is then -to a good approximation- 
replaced by 

^ 1 cosh Y± ^^g^ 

2 cosh^ ^7 [1 + sinh^ Y± tanh^ r]]^/^ 

where Y± is the transverse rapidity, Y± = 0.51og[(l + v)/{l — v)], and v is 
the tranverse velocity of the cluster. 

We have calculated the width of the balance function with ()16|) approx- 
imated by a Gaussian^ of the same width as that of (fT^ . The results are 
shown in Fig.l where < \6\ > is plotted versus v, for A = 3.5. The measured 
values for most central events (as reported in 2j) are also indicated. One 
sees that the calculated width decreases with increasing transverse velocity 
of clusters. One also sees that to obtain quantitative agreement with data 
the transverse velocity must approach 0.8c. Such a large value seems difficult 
to reconcile with other estimates of the transverse velocity |10j . 

6. The pion cluster model discussed so far ignores entirely the parton 
structure of the final system of hadrons. One may therefore not be surprized 
that it fails to describe the data from central heavy ion collisions. This 
argument suggests to try a model in which the parton structure is built in 
from the beginning. To this end we investigated coalescence model [H] which 
we generalized to include the correlations inside the system. 

To introduce correlations we assume that -just before hadronization- the 
QGP forms the weakly correlated neutral clusters. The clusters decay into 
quarks'^, antiquarks and gluons. One cluster provides one^ qq pair (either uu 

®It was shown in 9^ that this is a good approximation. 
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Figure 1: Width of balance function versus velocity of transverse flow. 



or dd) and any number of gluons. In the final step quarks and antiquarks 
coalesce into observed hadrons. The remaining gluons form again neutral 
clusters and the process continues. 

Thus the model we consider is basically the well-known coalescence model 
jH] supplemented by a prescription for correlations. Since the coalescence 
model was rather succesful in description of single particle spectra in central 
collisions of heavy ions [HIE]; it seems worthwhile to investigate its extention 
to correlation phenomena (see also |I2])- Admittedly, the proposed extension 
is very simple - perhaps even simplistic. It contains, however, all ingredients 
necessary to formulate and study the width of balance functions which is of 
interest in this paper. Therefore we do not find necessary at the moment to 
formulate and discuss a more general and/or detailed approach. 

7. To evaluate the balance function we need the distribution of pairs of 
charged pions, same charge as well as opposite charge. The pairs of same 
charge can be constructed by coalescence of the decay products of four clus- 
ters (two U-clusters and two D-clusters)''. The distribution of the pairs of 

^To shorten the wording, we caU by U-cluster the one decaying into uu and by D-cluster 
the one decaying into dd. Their distributions and decay properties are identical. 
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opposite charge consists of two terms: one identical to the distribution of 
same charge pairs and another one, arizing from coalescence of the decay 
products of two clusters (one U-cluster and one D-cluster). Thus the contri- 
butions involving four clusters exactly cancel and we only have to consider 
the distribution of pions of opposite charge which result from coalescence of 
decay products of one U-cluster and one D-cluster. 
This distribution can be expressed as 

Piv+,V-) = I dYudYDPciYc^Y) 

dT]udT]ufq{T]u - Yu) fq{r]u - Yu)hg{[r]u + 7]u]/2 - Yu) 

dVddr]afgivd - YD)fqiVd - YD)hg{[r]d + r]^]/2 - Yd) 

- ijiu + m)l'A^[v^ - ivd + Vu)/2Mvu - vMvd - Vu] (i7) 

where fq and hq are responsible for the distribution of quarks in decay of 
either U or D cluster, while $ summarizes the properties of the coalescence 
process. Finally, paiXc, ^y) denotes the joint distribution of U and D clus- 
ters with average rapidity Yq, where 

>c = ^^^; Ay = Yu-Yn (18) 
To simplify the discussion, in the following we shall assume that pc factorizes: 

Pg{Yc,Ay)=Pc{Yc)p{Ay) (19) 
Taking advantage of the delta functions we can rewrite ()17p as 



piv+,V-) = J dYcdAYPGiYcAy) J 



duj^du- 



U (.. + ^ - Yc) U (.- + ^ - Yc) K (. + ^ - Yc 
U (V - \ - Yc) U (.. - ^ - Kc) K \ - ^ - Ya 

$ (m+ + Ay) $ (U_ + Ay) (20) 

where z = (//+ + ri-)/2. 

To proceed, we again consider Gaussians 

fq{x) = ^e-"'/^'; hq{x) = Jl + ay2h^exp[-xyh^] 

$(x) = ^e-^lp\ (21) 



With this Ansatz, the integrals over du^du^ can be performed. The result 

is 

p{v+, V-) = Ce-'"!'" J dYcPc{Yc) exp 

where S = t]^ — rj^ and C is a constant, irreleveant for further discussion. 

The formula (j22p can be now introduced into (0) and thus the balance 
function can be calculated. In the limit of very large acceptance we obtain 



A[z-Yc 



c 



■{1 + ^/2^) 



(22) 



5s(5;A)A^oo = ^e-'^'/^' (23) 

One sees that -in this limit- the width of the balance function depends on 
one parameter which -to a large degree- determines also the distribution in 
decay of a cluster into free quark and antiquark. Indeed, using (fTTj) and (EU, 
one can show that the decay distribution in the rest frame of the cluster is 
given by 



PpiVu) = 7^ exp 



rf2 



(24) 



where 



2l + cV4/l2 



d' = c' ' (25) 

8. It seems natural to assume that -in their rest frame- clusters decay 
isotropically. This means that their decay distribution is given by (jl5j) . the 
same as for the clusters of pions considered before. It follows that -for an 
uncorrelated decay {cjh ^ 0)- the parameters / in (fTIH) and c in (PTj) are 
identical. Comparing (fT^ and ()23|) we thus conclude that in the coalescence 
model the width of the balance function is expected to be by factor a/2 smaller 
than that obtained for pion clusters. The reason is clear: the dispersion of 
the pion rapidity is reduced by precisely this factor when the pion is formed 
by random coalescence of a quark and an antiquark. 

Repeating the argument of the section 5, we thus conclude that -ignoring 
for the moment the corrections for finite acceptance and effects of transverse 
flow- the width of the balance function is expected to lie between .69 and .98 
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(the lower limit is obtained for a = {c/hY = 0, i.e. when the decay products 
of a cluster are uncorrelated) . 

To compare this result with the data we have to estimate the corrections. 
To this end we take the Gaussian Ansatz for pc'- 

PciYc) = -T^e-'^S/^' (26) 
which allows to evaluate explicitely the integrals in (j22|) . We obtain 



where 



c 



+ 4^2(1 + cV2/i2 



Using ()27|1 one can now follow the argument of section 5 and evaluate the 
width of the balance function, taking into account the finite acceptance and 
the transverse flow. In Figure 1 the width of the balance function evaluated 
from (P7|) is plotted versus f for A = 3.5 and two values of the parameter a = 
c^/h"^. One sees that these effects reduce substantially the calculated width. 
The value found in [2] for central collisions is reproduced with transverse 
velocity below 0.5, consistent with other estimates of the transverse flow 

[mi. 

One also sees from the Fig. 1 that in the coalescence model the calculated 
width is smaller than the value 0.65 found in j2j for peripheral collisions. This 
is not surprizing: in peripheral collisions a substantial part of the particle 
production should resemble the elementary nucleon-nucleon collisions which 
are not expected to follow the coalescence mechanism yjjj and are charac- 
terized by a significantly larger width of the balance function PJEj. As seen 
from Fig.l, the width of the balance function calculated from the pion cluster 
model (adequate for nucleon-nucleon collisions) is indeed close to 0.65. 

9. In conclusion, we have shown that the coalescence mechanism implies 
a substantial reduction of the pseudorapidity width of the balance function. 
This allows to explain the small width observed for central collisions of heavy 
ions [2:, provided the corrections due to the finite acceptance region and to 
the transverse flow are taken into account. This result supports the coales- 
cence mechanism as the final stage of the process of hadronization. 
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